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Abstract 

We consider the Dirac equation in with constant coefficients and study the 
distribution fit of the random solution at time t G IR . It is assumed that the ini- 
tial measure /^o has zero mean, a translation-invariant covariance, and finite mean 
charge density. We also assume that /^o satisfies a mixing condition of Rosenblatt- 
or Ibragimov-Linnik-type. The main result is the convergence of to a Gaussian 
measure as t — >■ c« . The proof uses the study of long time asymptotics of the solution 
and S.N. Bernstein's "room-corridor" method. 
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1 Introduction 



This paper can be regarded as a continuation of our papers [l]-[5] concerning the analysis 
of long-time convergence to an equilibrium distribution for hyperbolic partial differential 
equations and harmonic crystals. Here we develop the analysis for the Dirac equation 

t) = [—a ■ V — iPm] t), x E JR^, ^ 

where V = 82, d^) , — d/dx^ , m > , a — (cci, 0:2, 0:3) , and (3 are 4x4 Dirac 
matrices (see (2.12), (2.13)). The solution ^{x,t) E for {x,t) e IR^ . 

It is assumed that the initial data ipoi^) ^'I's given by a random element of the function 
space H = Hi^^ClR^) of states with finite local energy, see Definition 2.1 below. The distri- 
bution of ■00 is a zero-mean probability measure /xq satisfying some additional assumptions, 
see Conditions S1-S3 below. Denote by Ht , i G IR , the measure on H giving the distribu- 
tion of the random solution ilj{t) of problem (1.1). We identify the complex and real spaces 
(C^ = IR^ , and stands for the tensor product of real vectors. The correlation functions of 
the initial measure are supposed to be translation-invariant, 

Qo{x,y) -.^ E{7j;o{x)^^o{y)) ^qo{x-y), x,yeJR^. (1.2) 

We also assume that the initial mean charge density is finite, 

eo := £^|V'o(2;)|^ = trgo(O) < 00, xeJR^. (1.3) 

Finally, assume that the measure /Iq satisfies a mixing condition of a Rosenblatt- or Ibragimov- 
Linnik type, which means that 

ipo{x) and ipo{y) are asymptotically independent as \x — y\ ^ 00. (1.4) 

Our main result gives the (weak) convergence of /it to a limit measure /^oo , 

Ht^l^oo: t^oo, (1.5) 

which is an equilibrium Gaussian measure on 7-^ . A similar convergence holds for t ^ —00 
because our system is time-reversible. Explicit formulas (2.17) for the correlation functions 
of //(X, are given. 

To prove the convergence (1.5) we follow the strategy of [l]-[5]. There are three steps. 

I. The family of measures /if , t > , is weakly compact in an appropriate Prechet space. 

II. The correlation functions converge to a limit, 

Qtix,y) = J ip{x) <S)ijj{y) HtW ^ Qoc{x,y), t ^ 00. (1.6) 

III. The characteristic functionals converge to a Gaussian functional, 

fiM) := J e*<^'^>/.t(d^) ^ exp{-^Qoo(0, 0)}, t ^ 00. (1.7) 

Here (•, •) stands for a real scalar product in , for a quadratic form with the integral 
kernel Qoo{x,y) , and for an arbitrary element of the dual space. 
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Property I follows from the Prokhorov Compactness Theorem by a method used in [13]. 
Namely, we first establish a uniform bound for the mean local charge with respect to the 
measure //j , t > 0. Then the Prokhorov condition follows from the Sobolev embedding 
theorem by Chebyshev's inequality. Property II is derived from an analysis of oscillatory 
integrals arising in the Fourier transform. However, the Fourier transform by itself is insuffi- 
cient to prove Property III. We derive it by using an explicit representation of the solution in 
the coordinate space with the help of the Bernstein' "room-corridor" technique by a method 
of [l]-[5]. The method gives a representation of the solution as a sum of weakly dependent 
random variables. Then (1.5) follows from the Ibragimov-Linnik central limit theorem under 
a Lindcbcrg-type condition. Wc sketch the proofs by using the technique of [1]. 

The paper is organized as follows. The main result is stated in Section 2. The compact- 
ness (Property I) is estabfished in Section 3, the convergence (1.6) in Section 4, and the 
convergence (1.7) in Sections 5. 

2 Main results 

Let us describe our results more precisely. 
2.1 Notation 

We assume that the initial date ipo in (1.1) is complex- valued vector function belonging to 

the phase space Ti. . 

Definition 2.1 Denote by H = Hll^^(]R'^, (D"^) the Frechet space of complex-valued functions 
il!{x) , endowed with local energy seminorms 

\mlR= I mx)\'dx<oc, yR>0. (2.1) 

|a;|<i? 

Proposition 2.2 (i) For any ipo & H there exists a unique solution ip{-,t) e C(IR, H) to 

Cauchy problem (1.1). 

(ii) For any t G IR, the operator U{t) : ipQ i— > ip{-,t) is continuous in Ti . 

Proposition 2.2 follows from [10, Thms. V.3.1, V.3.2]) because the speed of propagation 
for Eq. (1.1) is finite. 

Let us choose a function C,{x) e C^(IR^)such that C(0) ^ 0. Denote by //f„^(IR^), 
s e H, the local Sobolev spaces, i.e., the Frechet spaces of distributions u e D'(IR^) with 
the finite seminorms 

\\u\1,r:= ||A^(C(^/^)«)||l^(]r3), (2.2) 
where A^^; := F^}^^{{kYv{k)) , {k) := ^J\i^ -\- 1 , and i) :— Fv is the Fourier transform of 
a tempered distribution v . For (f) e C^(IR^) write F(f){k) = j e^'''''(f){x)dx. Note that the 
space ifjo(,(IR^) for s = agrees with Definition 2.1. 
Definition 2.3 For seM, write W = Hf^^(5R^). 

Using the standard technique of pseudodifferential operators and Sobolev's embedding 
theorem (see, e.g., [8]), one can prove that H — TiP C TiT^ for every £ > 0, and the 
embedding is compact. 
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2.2 Random solution. Convergence to equilibrium 

Let (n, S,P) be a probability space with expectation E and let BiTi.) be the Borel cr- 
algebra of TC. Assume that ipo = ipoi,'-^,-) iii (1-1) is a measurable random function with 
values in (H, BiTi)) . In other words, (a;, x) i— > ^) is a measurable mapping $1 xIR^ — > 
with respect to the (completed) cr-algebras E x B(IR^) and B(C^) . Then, by virtue 
of Proposition 2.2, ilj{t) = U{t)il)Q is again a measurable random function with values in 
(7i, B{T-C)) . Denote by fioldipo) the Borel probability measure on Ti giving the distribution 
of ipo . Without loss of generality, we can assume that (f2, E, P) = {7i, B{7i), Ho) and 
11)0(01, x) = uj{x) for iJ,o{duj) X da; -almost all points {uj,x) e ?i x IR^ . 

Definition 2.4 Let Ht be the probability measure on Ti giving the distribution of Y{t) , 

lit{B) = i^o{U{-t)B), VP e B{n), t e JR. (2.3) 

Our main objective is to derive the weak convergence of the measures /if in the Prechet 
space for each £ > , 

fJ-t — r as t^ 00, (2.4) 
where //oo is some Borel probability measure on the space H . This means the convergence 

J fWMdi^) ^ / /(^)/ioo(#), t ^ 00, (2.5) 

for any bounded continuous functional /('?/') on TC~^ . 

Set TZip = (Re-?/', Im-?/') = {Re-i/'i, . . . , Re'?/'4, Im-j/'i, . . . , Im-?/'4} for ip = {ipi, . . . e 
, and denote by Wip j th component of the vector TZijj , j = 1, 8 . The brackets (•, •) 

mean the inner product in the real Hilbert spaces = L^(IR^) , in (8) IR^ , or in some 

their extensions. For iIj{x), (f){x) e L^(IR^, (D^) , write 

:= (7^V',7^0) = ^(7^^v^,7^^0). (2.6) 

Definition 2.5 The correlation functions of the measure /it are defined by 

Q"/{x,y) = E{Wip{x)Wil){y)) for almost all x,y e 1R^ i,j = 1,...,8, (2.7) 
provided that the expectations in the RHS are finite. 

Denote by D the space of complex-valued functions in C^(1R'^) and write T> := [P]^ . 
For a Borel probability measure n on 7i , denote by fi the characteristic functional (the 
Fourier transform) 

/i(0) = yexp(i(^,(/)))/x(#), <peV (see (2.6)). 

A measure /i is said to be Gaussian (with zero expectation) if its characteristic functional 
is of the form 

/i(0)=exp{-^Q(0,0)}, 0eP, 

where Q is a real nonnegative quadratic form on P . A measure /i is said to be translation- 
invariant if 

fi{ThB) ^ n{B), BeB{n), heJR\ 
where Thijj{x) = ijj{x - h) , x eJR^ . 
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2.3 Mixing condition 



Let 0(r) be the set of all pairs of open bounded subsets A.Bc IR'^ at the distance 
dist {A, B) > r , and let (j{A) be the a -algebra in 7i generated by the linear functionals 
■01-^ , where (f) eT> with supp (p G A. Define the Ibragimov-Linnik mixing coefficient 

of a probability measure //q on H by formula (cf. [9, Def. 17.2.2]) 

ip{r) = sup sup — . (2.8) 

(AB)GO(r) AeaiA),B ea{B) ^^^^> 

MB) > 

Definition 2.6 We say that the measure /iq satisfies the strong, uniform Ibragimov-Linnik 
mixing condition if 

(p{r) ^0 as r — > CX3. (2.9) 
We specify the rate of decay of ip below (see Condition S3). 

2.4 Main assumptions and results 

We assume that the measure no has the following properties SO— S3: 

50. Ho has zero expectation value, EiPq{x) = 0, x e IR^ . 

51. //o has translation-invariant correlation functions, 

Q^{x,y) = E(n'ip{x)Wip{y)) = qo^{x-y) for almost all x,y e 1R^ i,j = 1, ...,8. (2.10) 

52. fiQ has a finite mean charge density, i.e., Eq. (1.3) holds. 

53. Ho satisfies the strong uniform Ibragimov-Linnik mixing condition, with 

/•oo 

/ r^(p^/'^(r)dr < oo. (2.11) 
Jo 

The standard form of the Dirac matrices ctfe and /3 (in 2 x 2 blocks) is 
where / denotes the unit matrix and 

Introduce the following 8x8 real valued matrices (in 4x4 blocks) 



Note that by (2.12) and (2.13) we have 

ia2 



ia2\ , . / 1 
ta2 = . , where ta2 



Moreover, Al = Ak , k = 1,2,3, Aq = -Aq • Write 



A = (Ai, A2, A3), P(V) = A • V + mAo. 



(2.15) 



For almost all x,y E IR^ , introduce the matrix- valued function 



Q^{x,y) = {QU^,y)),.^,_, = (?^(^ - y)) 



(2.16) 



Here 



qoo{k) = + \v{k)P{-ik)%{k)P'^{ik), 



(2.17) 



where V{k) 
the measure 



= 1/(A;^ + m^) , and qo{k) is the Fourier transform of the correlation matrix of 
jiQ (see (2.10)). Since P'^{ik) = —P{—ik) , we have, formally. 




where V{z) = e~"^l^l/(47r|z|) is the fundamental solution for the operator —A + , and 
* stands fors the convolution of distributions. We show below that % e L'^ = L^(IR^) (cf 
(4.7)). Hence, qoo{k) G by (2.17), and the convolution in (2.18) also belongs to . 

Denote by Qoo a real quadratic form on defined by 



The form Qoo is continuous in because ^00 (^) is bounded by Corollary 4.3. 

Theorem A. Let m > , and let S0-S3 hold. Then 
(i) the convergence in (2.4) holds for any e > . 

(a) The limit measure fi^o is a Gaussian equilibrium measure on Ti. . 
(Hi) The characteristic functional of /loo is of the form 



Theorem A can be derived from Propositions 2.7 and 2.8 given below by using the same 
arguments as in [13, Theorem XII. 5. 2]. 

Proposition 2.7 The family of measures G IR} is weakly compact in the space Ti"^ 

for any £ > . 

Proposition 2.8 For any cj) e V , 



Q^{(t>,(t>) = {Q^{x,y),n(t>{x)®n(t>{y))= E / Q'i{x,y)1V(t>{x)W(t>{y) dxdy 





Propositions 2.7 and 2.8 are proved in Sections 3 and 4-5, respectively. 
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2.5 Remark on various mixing conditions for the initial measure 

Wc use the strong uniform Ibragimov-Linnik mixing condition for the simphcity of our 
presentation. The uniform Rosenblatt mixing condition [12] with a higher degree > 2 in 
the bound (1.3) is also sufficient. In this case we assume that there exists a S , S > , such 
that sup E\ipo{x)\'^'^^ < oo. Then condition (2.11) requires the following modification: 

roo 

/ raP (r)dr < OO, p = min((5/(2 + 5), 1/2), 
Jo 

where a{r) is the Rosenblatt mixing coefficient defined as in (2.8), but without the denomi- 
nator Ho{B) . The statements of Theorem A and their proofs remain essentially unchanged. 

3 Compactness of measures 

3.1 Fundamental solution of the Dirac operator 

One can easily check that au and /3 are Hermitian symmetric matrices satisfying the anti- 
commutation relations 

a* = ak, 13* = (3, 

akai + aittk = 26kil, ak(3 + (3ak = 0, 

( Ski is Kronecker's delta). Therefore, 

(dt + a-V + i(3m) {dt-a-V - i(3m) = (d^ - A + m^)7. 

Then we can construct a fundamental solution £{x, t) of the Dirac operator, i.e., a solution 
of the equation 

{dt + a-V + iPm)£{x,t) = S{x,t)I, £{x,t)^0 for t < 0, 

of the form 

£{x,t) = {dt- a-V - i(3m) E{x,t), (3.1) 

where E{x,t) = Et{x) is a fundamental solution for the Klein-Gordon operator {d^ — A -|- 
m^) , and E vanishes for i < . 

Remark 3.1 The function Et{x) is given by 

Et{x)^F,-^f^, u;^u;{k)^^\k\^ + m^. (3.2) 

Then, by the Paley-Wiener Theorem (see, e.g., [6, Theorem II. 2. 5.1]), the function of Et{-) 
is supported by the ball \x\ < t . 

Denote by U{t), t G IR, the dynamical group for problem (1.1). Then U{t) is a 
convolution operator given by 

V'(x, t) = [/(i)V'o = £{■, t)*ijjo^{dt-a-V - i(3m) Et{-) * V'o- (3.3) 

The convolution exists because the distribution £{■, t) is compactly supported by (3.1) and 
by Remark 3.1. 
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3.2 Local estimates 

Proposition 3.2 For every i/jq e H and R> , 

\\Umo\\o,R<C\Mo,R+t, teJR, (3.4) 
where C < oo does not depend on R and t . 

Proof. In the Fourier transform, the solution il}{x,t) of the Cauchy problem (1.1) reads as 

■0(/c, t) — S{k, t)ijjo{k) — cosut — (a ■ (—ik) + i(3m) '0o(^) 
by (3.1) and (3.2). Then, for V'o e , 

M-Ml^ = M-Ml^ < C\\M-)\\l^ = C'||V'o(-)IU- (3.5) 

Let us consider ipQ & Ti . Introduce the function ipo{x) equal to ipo{x) for < R + t 
and to otherwise. Denote by ip{x,t) (by ip*{x,t)) the solution of the Cauchy problem 
(1.1) with the initial data i/joix) ( i/^oix) , respectively). Note that ii}{x,t) = '4>*{x,t) for 
\x\ < R. Therefore, relation (3.5) implies 

m-Mn = wn-MR < c\m-)\\L^ = ^^ii^o(-)iu+*- ° 

3.3 Proof of compactness 

Proposition 2.7 follows from the estimate (3.3) below by using the Prokhorov Theorem [13, 
Lemma II.3.1], as in the proof of [13, Thm. XII. 5. 2]. 

Proposition 3.3 Let the conditions of Theorem A hold. Then, for any positive R, there 
exists a constant C{R) > such that 

snp E\\U(t)'ijjo\\lR < C{R) < oo. (3.6) 

t>o 

Proof. Let us write 

et{x) -.^ E\'il){x,t)\^, xelR^ (3.7) 

The mathematical expectation is finite for almost every x by (3.4) and by the Fubini the- 
orem. Moreover, et{x) = Ct for almost every x E by Condition SI. Hence, it follows 
from the Fubini theorem, (3.4) and Condition S2 that 

E\\U{t)^o\\l,R = et\BR\ < C£;||V'o||o,R+|t| = Ceo\BR+\t\l t e IR. (3.8) 

Here Br is the ball < R in IR^ , and \Br\ is the volume of this ball. As — > cxo , we 
see from (3.8) that Ct < Ccq . Thus, 

E\\U{t)MlR = ^ <^^o\Br\ < oo. □ 
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4 Convergence of correlation functions 

We prove the convergence of the correlation functions for the measures Ht ■ This imphes 
Proposition 2.8 in the case of Gaussian measures /iq ■ It follows from condition SI that 

Q?{x,y) = qi:^{x-y), x,yeJR\ (4.1) 

for i,j = 1,...,8. 

Proposition 4.1 The correlation functions qt{z), i,j — 1,...,8; converge for any z e 

IR^ 

qn^)^qZiz), t^oo, (4.2) 
where the functions 5^(2;) are defined in (2.17). 

Proof. Using the notation (2.14) and (2.15), by (3.3) we obtain 

nil^i^x, t) = (dt - P{V))Et * Ui/Jo. 

Then, by (3.2) and (2.15), the Fourier transform of the solution to Cauchy problem (1.1) 
becomes 

'R^{k,t) ^gt{k)'R^o{k), where Gtik) -.^ cos out- P{-ik)^^^^. (4.3) 
The translation invariance (2.10) implies that 

E^R^oik) (8) TZi^oik')) = F^^k,y-.k' Qoi^ - y) = (27rf 5(k + k')qo(k). (4.4) 
Further, (4.3) gives 

E{1Zij{k, t) 1^{k', t)) = {27rf 5{k + k')Gt{k)qQ{k)G*{k). (4.5) 
Therefore, by the inverse Fourier transform we obtain 
qt{x-y) = E{n^l;{x,t)®n^l:{y,t))^F-}^^^^_Jt{k)Uk)Q;{.k) 

= (27r)-3 /e-^'=(^-^)(cosc^i-P(-ifc)^)go(A;)(cosa;t-P^(ifc)^)dA; 
= (2.)-' / e--"-»> - 5i^(j„(*.)f.-(«) + P[-rkmk)) 

+ *• (4-6) 

To prove (4.2), it remains to show that the oscillatory integrals in (4.6) converge to zero. 
Let us first analyze the entries of the matrix i '^i3 — ^1 ■■■1'^ ■ 

Lemma 4.2 Let the assumptions of Theorem A hold. Then q^ e L^(IR^) for any i,j . 
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Proof. Let us first prove that 

ci^{z) e Lf(IR3), p>l, i,j = l,...,8. (4.7) 
Conditions SO, S2 and S3 imply by (cf. [9, Lemma 17.2.3]) that 

\ci^{z)\<Ceo^^'Wz\), ^elR^ z,j = l,...,8. (4.8) 
The mixing coefficient (/? is bounded, and hence relations (4.8) and (2.11) imply (4.7), 

oo 

j \(i^{z)\Pdz < Cel J </^'(l-2|) dz<CiJ r^^^/^r) dr < oo. 

]R3 R3 

By Bohner's theorem, is a nonnegative matrix-valued measure on IR^ , and condition S2 
implies that the total measure go(IR^) is finite. On the other hand, relation (4.7) for p — 2 
gives qi)^ G L'^{'iR^) . Hence, g^^' e L^(IR^) . □ 

Let us apply this lemma to the oscillatory integrals entering (4.6). The convergence (4.2) 
follows from (4.6) by the Lebesgue-Riemann theorem. This completes the proof of Proposi- 
tion 4.1. □ 

Relation (4.7) for p = 1 implies now that qo{k) is bounded. Hence, the explicit formula 
(2.17) implies the following assertion. 

Corollary 4.3 All matrix elements q]^{k) , i, j = 1, . . . ,8 , are bounded. 

5 Convergence of characteristic functionals 

To prove Proposition 2.8 for the general case of a non-Gaussian measure Hq , we develop a 
version of Bernstein's "room - corridor" method of [l]-[4]: (i) we use an integral represen- 
tation for the solutions of (1.1), (ii) divide the region of the integration into "rooms" and 
"corridors" and (iii) evaluate their contribution. As the result, the value {U{t)iljQ,(j)) for 
(p E V is represented as the sum of weakly dependent random variables. Then we apply 
Bernstein's "rooms-corridor" method and the Lindeberg central limit theorem. 

(i) We first evaluate the inner product {U{t)ipQ, 0) in (2.19) by using duality arguments. 
For t G IR , introduce a "formal adjoint" operators U'{t) from the space X> to a suitable 
space of distributions. For example, 

{^,U'{t)<l>)^{U{t)^,<l>), <f>ev, i;en. (5.1) 
Write 0(-,t) = U'{t)(j) . Then (5.1) can be represented as 

(V^(i),0) = (^o,0(-,i)), ielR. (5.2) 

The adjoint groups admit a convenient description (see Lemma 5.1 for the group U'{t) ). 
Lemma 5.1 For (f) &V , the function U'{t)(f) = (j){x,t) is the solution of 

4>{x,t) = {a -V + il3m)(f){x,t), 4>{x,0) ^ (f){x). (5.3) 



9 



Proof. Differentiating (5.1) with respect to t for ip,(j) eV , we obtain 

{^l;,U'{t)<l>)^{U(t)^,<l>). (5.4) 

The group U{t) has the generator A — —a • V — i/3m. Therefore, the generator of U'{t) is 
the conjugate operator 

A' ^a-V + i(3m. (5.5) 
Hence, relation (5.3) holds indeed with = A!(t). □ 

Remcirk 5.2 Comparing (5.3) and (1.1), we see that (f){x,t) = U'{t)(f) can he represented 
as a convolution (cf (3.3)), namely, 

ci>{;t) = nt*cP, nt:=idt + a-V + imP)Et. (5.6) 

(ii) Introduce a "room-corridor" partition of IR'^ . For a given t > , choose dt > I and 

Pt > such that pt ~ t^~^ with some 6 e (0, 1) and df ~ t/lnt , as t — >• oo. Set ht = dt + pt 
and 

a^^jht, V^a^ + dt, jeTL. (5.7) 

We refer to the slabs El = {x e IR^ : d' < < V} as "rooms" and to C{ ^ {x E 
JR^ : b' < x^ < %+i} as '^corridors" . Here x = {x^,x'^,x^) , the symbol dt stands for the 
width of a room, and for that of a corridor. 

Denote by Xr the indicator of the interval [0, dt] and by Xc that of [dt, ht] , which means 
that 

"^{Xris — jh) + Xc{s — jh)) ~ 1 for (almost all) s e M. 
The following decomposition holds: 

(^0, </>(-, t)) = E((^o,X^.</'(-,t)) + (^o,X^c'/'(-,^))), (5-8) 

where := Xr(2^^ ~ jh) and := Xc(2^^ ~ jh) ■ Consider the random variables rl and 
given by 

W = (V'o,X^.0(-,i)), 4 = (V'o,X^c0(-,i)>, (5-9) 
Then (5.8) and (5.2) imply 

(C/(i)V'o,0) = E(^t+c^')- (5-10) 

i6Z 

The series in (5.10) is in fact a finite sum. Indeed, for the support of we have 

supp (p G Br for some r > 0. 

Then, by the convolution representation (5.6), the support of the function (/){•, t) at i > 
is a subset of an "inflated future cone" 

supp(/> C {(.T,t) G IR^ X IR+ : |a;|<t + r}, (5.11) 

whereas 7lt{x) is supported by the 'future cone' |a;| < t. The latter fact follows from (5.6) 
and from Remark 3.1. Finally, it follows from (5.9) that 

rl — Ct for jht + t < —r and for jht — t>r. (5-12) 
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Therefore, the series (5.10) becomes a sum, 



Nt f 



{Uit)^l;o,<i>) = J2i'i + ci), Nt^-. (5.13) 



Lemma 5.3 Let Conditions SO— S3 hold. Then the following hounds hold for t > 1 : 

E\ri\'<C{cj>) dt/t, E\4\'<C{cj>) p,/t, jeZ. (5.14) 

Proof. We discuss the first bound in (5.14) only, because the other can be proved in a 
similar way. Rewrite the LHS of (5.14) as the integral of correlation functions. We obtain 

E\ri\^ = {xi{x3)xi{y3)qo{x - y),(j){x,t) ^ (j){y,t)). (5.15) 

The following uniform bound holds (cf. [11, Thm. XI. 17 (b)]): 

sup |(/)(x,i)| = 0(r^/2), t^oo. (5.16) 

In fact, (5.6) and (3.2) imply that the function (f){x,t) can be represented written as the 



sum 



cp{x,t) = E / e-^('=^^^*^a±(cu)0(A;) dk, (5.17) 



where a (u) is a matrix whose entries are hnear functions of a; or 1/uj . Let us prove the 
asymptotics (5.16) along each ray x — vt + Xq with \v\ < 1 . The asymptotic relation thus 
obtained must hold uniformly in x e IR^ by (5.11). By (5.17) we have 

cf>{vt + xo,t)^Yl J e-^('="^'")*-^'=^°a±(a;)0(A;) dk. (5.18) 

This is a sum of oscillatory integrals with the phase functions (f)±{k) = kv ± uj{k) . Each 
function has two stationary points which are solutions of the equation v = ^Va;(A;) if 
It"! < 1 , and has none if 1^1 > 1 . The phase functions are nondegenerate, i.e.. 

Finally, 0(/c) is smooth and rapidly decays at infinity. Therefore, (f){vt + xo,t) = 0{t^^^^) 
according to the standard method of stationary phase, see [7] . 
According to (5.11) and (5.16), it follows from (5.15) that 

E\ri\^<Ct-^ I xi{x')\\qo{x-y)\\dxdy = Ct-^ J xi{x')dx J\\qoiz)\\dz, (5.20) 

\x\<t+r \x\<t+r r3 

where ||?o(-z)|| stands for the norm of the matrix (^^{zfj . Therefore, relation (5.14) follows 
for \\qo{-)\\ e L\-\R') by (4.7). □ 

Hence, the rest of the proof of Proposition 2.8 is the same as that in the case of the 
Klein-Gordon equation, [1, p. 20-25]. The proof of Theorem A is complete. 
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